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- , Abstract: The effect of an external (electro) magnetic field on the finite temperature transition of 
O . QCD is studied. We generate configurations at various values of the quantized magnetic flux with Nf = 



2 + 1 flavors of stout smeared staggered quarks, with physical masses. Thermodynamic observables 
including the chiral condensate and susceptibility, and the strange quark number susceptibility are 
measured as functions of the field strength. We perform the renormalization of the studied observables 
and extrapolate the results to the continuum limit using Nf = 6,8 and 10 lattices. We also check 
for finite volume effects using various lattice volumes. We find from all of our observables that the 
transition temperature Tc significantly decreases with increasing magnetic field. This is in confiict with 
various model calculations that predict an increasing Tc{B). From a finite volume scaling analysis we 
find that the analytic crossover that is present at i? = persists up to our largest magnetic fields 
eB ~ 1 GeV^, and that the transition strength increases mildly up to this eB ~ 1 GeV^. 
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1 Introduction 

The properties of QCD in strong magnetic^ fields are relevant for at least three important physical 
situations. First, cosmological models suggest that extremely strong magnetic fields [VeB ~ 2 GeV) 
could be produced during the electroweak phase transition of the early universe. This effect might also 
have an impact on subsequent strong interaction processes [1]. Second, large magnetic fields [V^B ~ 1 
MeV) are present in the interior of dense neutron stars called magnetars [2]. Finally, in a noncentral 
heavy ion collision the spectators - being two beams of positive charges moving in opposite directions 
- also create an intense magnetic field which, depending on the centrality and the beam momentum, 
reaches up to V eB ~ 0.1 GeV for RHIC and V eB ~ 0.5 GeV for the LHC [3]. This magnetic field is 
^Throughout the paper 'magnetic' refers to electromagnetic i.e. not chromomagnetic. 
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external since it is produced by the spectators, and though it has a very short hfetime (of the order 
of 1 fm/c), the magnetic 'impulse' coincides with the generation of the quark-gluon plasma and thus 
may have a significant effect on the properties of the transition. 

For noncentral heavy ion collisions, an exciting consequence of the interplay between the strong 
magnetic field and the nontrivial topological structure of the quark-gluon plasma is the so-called 
chiral magnetic effect [4, 5]. This effect creates an electric current of quarks (anti)parallel to the 
external magnetic field, which may result in a preferential emission of charged particles perpendicular 
to the reaction plane, leading to event-by-event CP-violation [6]. Recent measurements from the STAR 
experiment at RHIC [7, 8] and the ALICE experiment at the LHC [9] are in qualitative agreement 
with this picture, however, the interpretation of these results is still under discussion [10-12]. 

Because of this high phenomenological relevance, the effect of a finite magnetic field on the strong 
interactions has been studied extensively in the last years, both using model calculations and lattice 
simulations. In particular, the structure of the QCD phase diagram in the B — T plane has received 
increasing attention recently. Calculations have been carried out within various low energy effective 
models of QCD. In the linear sigma model coupled to quarks and the Polyakov loop it was observed that 
the transition temperature increases with B [13]. Furthermore, a splitting between the deconfinement 
and chiral transitions was predicted to take place for large external fields. The strength of the transition 
was also observed to increase, which eventually results in a first-order phase transition [14]. Similar 
conclusions with respect to the increase in Tc and in the strength of the transitions were drawn from 
studies of the Nambu-Jona-Lasinio model and extended versions thereof, like the EPNJL and PNJLg 
models [15, 16], see also [17], and the nonlocal PNJL model [18]. 

The presence of the external magnetic field was also shown to increase the transition temperature 
within other types of models like in the Sakai-Sugimoto model of large Nc gauge theories [19], in 
the Gross-Neveu model in lower dimensions [20, 21], in 2+1 dimensional QED when described by 
Schwinger-Dyson equations [22] and within the holographic approach [23]. However, the opposite 
effect of a decreasing deconfinement transition temperature was predicted using chiral perturbation 
theory for two quark flavors [24] . A decrease in Tc was also observed in the linear sigma model if the 
quark vacuum contributions are neglected [13] and in the Sakai-Sugimoto model with nonzero chemical 
potential [25]. We mention that lattice simulations indicate a reduction of the transition temperature 
of QCD in an external c/iromomagnetic field [26-28]. 

The phase diagram is in most cases predicted by studying chiral symmetry breaking, i.e. the 
behavior of the chiral condensate or the dynamical quark mass as a function of B. Most of the low 
energy models agree that chiral symmetry breaking is enhanced as the magnetic field B grows [29- 
31]; in particular the value of the chiral condensate was found to increase linearly with \B\ in leading 
order [32-34] . The condensate also increases with B - although with a quadratic leading order - within 
the AdS/CFT duality picture [35], and with B^^^ in holography [23]. On the other hand, it was also 
conjectured that the running of the strong coupling in the presence of magnetic fields may modify this 
magnetic catalysis, and even turn the effect around to make the dynamical mass decrease with B in 
some regions [36]. 

In recent lattice simulations with Nj = 2 flavors of staggered quarks [37] the chiral condensate 
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was observed to grow with the external field for any temperature T in the transition region. The size 
of this effect was however found to be different for different values of T, resulting in an increase in both 
the pseudocritical temperature Tc and the strength of the transition. Furthermore, according to the 
findings of [37], the relative change in Tc is of the order of one percent for several larger-than-physical 
pion masses. 

In this paper our aim is to perform a similar lattice study, but with improved gauge and smeared 
fermionic actions and with Nf = 2 + 1 flavors of quarks, at the physical pion mass, and extrapolate the 
results to the continuum limit. We include the magnetic field in the fermion determinant and study 
its effect dynamically to investigate how the strength and the pseudocritical temperature of the QCD 
transition change as the external magnetic field is switched on. We explore a wide temperature region 
around the zero-field pseudocritical temperature Tc{B = 0), for various values of the magnetic field, 
ranging from V eB ~ 100 MeV to \/ eB 1 GeV, i.e. covering the regions that are phenomenologically 
interesting for noncentral heavy ion collisions and for the evolution of the early universe. 

This paper is structured as follows: first the implementation of the magnetic field on the lattice is 
described. Then we define the observables of interest, including the chiral condensate, chiral suscepti- 
bility and strange quark number susceptibility, and discuss their renormalization at zero and nonzero 
B. After presenting the simulation setup and the details of the analysis we show our results for the 
transition temperature and the width of the transition. 

2 Magnetic field on the lattice 

Let us consider the case of a constant external magnetic field B = (0, 0,-B), that is pointing in the z 
direction. In the continuum such a magnetic field can be realized by, e.g., the following vector potential, 

= {A, At) = {0,Bx, 0,0). (2.1) 

Any other vector potential satisfying B = curl(A) corresponds to the same physical system, and is 
connected to the above choice by an appropriate U(l) gauge transformation. 

It is well known that in a finite box with periodic boundary conditions the magnetic flux cannot 
be arbitrary, but is quantized in terms of the area A of the system in the plane orthogonal to the 
external field [38, 39]. This leads to the quantization condition, 

qB-A = 2TrNb, Nb G (2.2) 

where q is the charge of the particle. In turn, on the lattice the area A is also quantized as A = NxNyO?, 
where a is the lattice spacing (we restrict the discussion to isotropic lattices) and Ny is the number of 
lattice points in the direction v. This implies that the lattice discretization also imposes an upper bound 
on the magnetic flux. To see this explicitly, let us write down how the continuum vector potential (2.1) 
can be represented by complex phases Uy{n) G U(l) that multiply the Uu{n) G SU(3) links of the 
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lattice, 



Uy{n) = e'" 

(2.3) 

Ux{n) = 1, Uxi^Nx- 1, 

where the sites are labeled by integers n = {nx,ny,nz,nt), with riy = . . . — 1. Constant magnetic 
background fields were first used in lattice studies of nucleon magnetic moments at zero tempera- 
ture [40-42]. At finite temperature this approach was first realized in [37, 43]. 

On the lattice the quantization condition (2.2) thus takes the same form as in the continuum with 
the area discretized in terms of the lattice spacing, 

In this formulation we have periodic boundary conditions in all spatial directions and the magnetic 
flux going through any plaquette in the x — y plane is constant. Furthermore, this implementation of 
the magnetic field is periodic in the fiux quantum A*";, with a period of NxNy, so its value is effectively 
constrained to < A^f, < NxNy. This prescription is discussed in more detail in appendix A. We note 
that the periodicity of the field in Nf, implies a stronger constraint for the fiux, 

< AT, < (2.5) 

where the correspondence between the implementation and the actual value of B is unambiguous. 
Namely, at larger values of A";, the periodicity is expected to introduce saturation effects, like it was 
observed in [44, 45]. The largest possible magnetic field is therefore qB"^^^ = 7r/2 • a~^. 

On the lattice the temperature of the system is given by the inverse temporal extension as T = 
(A^to)^^. It is therefore clear from equation (2.4) that the minimal value of the magnetic field is 
^^min _ rp2 . 2'KNf /NxNy. Thus, to increase the maximal field one has to decrease a, and to decrease 
the minimal magnetic field one has to increase NxNy. Taking these considerations into account, with 
reasonable lattice spacings and lattice extensions, the lattice magnetic field covers the region y/qB = 
0.1... 2 GeV. 

We remark that if there are particles with different charges in the system, then the quantization 
condition for B has to be fulfilled for the greatest common divisor - in our case this is the down quark 
charge, q = qd = —\e\/3 in equation (2.4) with e the charge of the electron. This will then determine the 
minimal field. Fortunately in nature the ratio of quark charges is a small natural number so that the up 
and down quarks can be studied together. We note furthermore that the above implementation of the 
magnetic field leads to no sign problem, in contrast to a finite chemical potential, or a (Minkowskian) 
electric field. 
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3 Observables at finite B 



Let us consider the staggered partition function with three flavors u,d and s. Each quark flavor has 
to be treated separately since the charges/masses are different: = —2qd = —2qs, and we assume 
viu = rrid 7^ ms- The partition function reads, after taking the fourth roots of the fermion determinants, 

VUe-P^o [det M(C/, g^, [det M{U, qa, m^, /i^)]'/^ [det M(C/, g,, m„ /i,)]^/^ , (3.1) 

where the fermion matrix is M(U, q, m, fi) = lJ){U, q, ^)+ml. (Here we do not address problems arising 
from the rooting trick [46]). The dependence on the chemical potential is only made explicit to define 
derivatives of the partition function with respect to see equation (3.3), and later we set all chemical 
potentials to zero. Since we are only concerned with a constant external field, the dynamics of the 
U(l) field introduced above does not have to be taken into account; in the gauge sector we only have 
the SU(3) kinetic term Sg with inverse gauge coupling 13 = . 

To study thermodynamics in a nonzero external field we analyze the chiral condensates and chiral 
susceptibilities for the light flavors f = u,d, 

- T dlogZ T dHogZ 
and the strange quark number susceptibility, 

T 1 dHogZ 

where we defined the spatial volume of the system as y = {NgaY with Ng = = Ny = N^- The 
condensate for a particular flavor will be denoted in the following by the first letter of the flavor name, 
e.g. uu. 

To take the continuum limit, the renormalization of these observables has to be carried out. The 
logarithm of the partition function log^ (i.e. the free energy) at i? = contains additive divergences 
of the forms a~^, nn?a~'^ and m^log(a) [47]. In section 4 we will show - based on the behavior of the 
beta function measured at zero temperature - that there are no additional i?-dependent divergences. 
Therefore the additive divergences of the observables derived from the free energy can be eliminated 
by subtracting the T = 0, i? = contribution. In the chiral quantities there are also multiplicative 
divergences caused by the derivative with respect to the quark mass. To eliminate this multiplicative 
divergence in the chiral condensate (susceptibility) , we multiply by the first (second) power of the bare 
quark mass [48]. Finally, to obtain a dimensionless combination we divide by the fourth power of the 
T = pion mass jn^, 

^PrfiB,T) = ruf [i^i^f{B,T) - i;^f{B = 0,T = 0)] ^, 

"'^ (3.4) 

x}iB, T) = m} [xf{B, T)-XfiB = 0,T = 0)] ^. 
Note that this procedure leads to a renormalized condensate that, for i? = 0, is zero at T = and 
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approaches a negative value as T is increased. 

Considering the strange quark number susceptibiHty, c| needs no renormahzation (neither at 
B = nor at i? 7^ 0) since it is connected to a conserved current. 

4 Renormalization at finite B 

We expect that a nonzero external magnetic field does not introduce new divergences in the free energy 
density, since the external field - just like a chemical potential - is coupled to the current ip'yui^ which is 
conserved. This expectation is also supported by the fact that in the presence of the external field there 
are no additional, divergent Feynman-diagrams due to the absence of internal photon lines. ^ Moreover, 
the vacuum energy was also calculated in the effective potential approach [49] and its divergent part 
was found to be independent of B. The absence (or presence) of i?-dependent divergences in the free 
energy density is closely related to the non-renormalization (or renormalization) of B itself. In fact, 
in a gauge invariant renormalization scheme the product eA^ needs no renormalization because of the 
U(l) Ward-Takahashi identity (see appendix B). In our case the magnetic field always appears in the 
combination eB and, therefore, we expect that it is not subject to renormalization"^. We now check 
these expectations numerically. 

From the point of view of renormalization theory it may be instructive to draw a parallel between 
the magnetic field and the quark mass, //the magnetic field were to induce new divergences (like the 
mass does) then eB itself (like m) would be subject to renormalization as {eBy = Z ■ (eB) with Z a 
corresponding renormalization constant. However, in the lattice approach the magnetic field - unlike 
the mass - fulfills a quantization condition, as in equation (2.4). Therefore the renormalization of eB 
can only amount to a shift in the lattice spacing a, such that a^eB — a^y^^f^^^leBY is satisfied. This 
implies that the lattice scale has to change if the renormalization of the magnetic field is nontrivial, 
i.e. if there are ei?-dependent divergences. This is expected since the lattice scale is determined by the 
beta-function of the theory which is given in terms of the renormalization scale-dependent, divergent 
Feynman-diagrams (see e.g. [50]). For the magnetic field however, due to the quantization condition, 
the only possible effect of such divergent diagrams is to alter the lattice scale. 

Therefore we propose to measure a physical quantity (/> at T = as a function of the magnetic 
flux Nf, for different lattice spacings. We take the lattice scale a(/3) and the line of constant physics 
(LCP) ?7i(/3) which are measured at Nf, = (see section 5), and assume that they are also valid at 
> 0. We use the scale to determine the magnetic field from the fiux according to equation (2.4) 
and the quantity <j) in physical units. Then we compare the running of (j){eB, a) with the magnetic 
field for different lattice spacings. If this quantity has a meaningful continuum limit, 

lim(/.(eS,a) =</.™^*(eS), (4.1) 

a— !>0 

^Consider e.g. the gluon self-energy diagram (with one quark loop) which is - in a gauge invariant regularization - 
logarithmically divergent. The coupling to the external magnetic field in the lowest order in B is given by two external 
photon legs attached to the quark loop. This diagram is clearly finite since it contains two extra quark propagators. 

^Note that for a dynamical U(l) theory, B would appear separately in the photon gauge action. This is, however, 
not the case for the present study. 
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i.e. if for small enough lattice spacings the dependence of (j){eB,a) on a is suppressed, then our 
assumption was valid and the lattice scale a(/3) and the LCP m(/3) are also correct for > 0. In view 
of the discussion in the previous paragraph, this suggests that there are no ei?-dependent divergences. 
In the opposite case the lattice scale does depend on N^, which in turn would indicate the presence of 
e-B-dependent divergences. 

For physical quantities we choose the charged pion mass (j) = and the Sommer parameter 

4> = tq. For their definition and measurement details see, e.g. [51]. For the charged pion mass we expect 
a strong dependence on the magnetic field, in the form [52], 



ml+iO) + \eB\, 



(4.2) 



which can be deduced from the dispersion relation for spin-1 mesons^. On the other hand, for the 
Sommer parameter, which is defined using the potential between static color charges, we expect the 
dependence on eB to be suppressed. 

In figure 1 the mass of the charged pion and the Sommer parameter are plotted as functions of 
the external field for various lattice spacings, at a = 0.29, 0.22, 0.15 and 0.12 fm (for the Sommer 
parameter the coarsest lattice spacing is not shown since here tq has large systematic errors). The 
lattice geometries and simulation parameters for these runs are tabulated in appendix C. At i? = we 
used the measurements presented in [51]. We observe a nice scaling with a for both quantities, in the 
region eB < 0.4 GeV^ where we have data for the three finest lattices. Results are consistent with a 
constant behavior for the Sommer parameter which indicates that the lattice spacing is not modified 
by the external field beyond our statistical accuracy. 
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Figure 1: The mass of the charged pion (left panel) and the Sommer parameter (right panel) as functions of 
the external magnetic field for different lattice spacings. Results for the pion mass are compared to the analytic 
prediction (see text). The good scaling of the lattice results and the independence of tq on eB indicate the 
absence of e_B-dependent divergences. 

Data for the charged pion mass are also as expected and agree with the analytic prediction (4.2) 
within 2 — 3%. For large Ni, ~ a^eB we see deviations from the continuum scaling only for the coarsest 
lattice which is most probably due to lattice artefacts stemming from the periodicity of the lattice 
magnetic field (2.3). Based on theoretical arguments (see appendix B) and on these observations we 



*We note that this expression ought to receive corrections for large B due to pair production. 
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conclude that it is safe to use the B = lattice scale and LCP at nonzero external fields, and to 
exclude the possibility of C-B-dependent divergences. 

5 Simulation details 

We use the tree-level improved Symanzik gauge action and stout smeared staggered fermions; details 
about the action can be found in [53]. We generate lattice configurations both at T = and T > with 
an exact RHMC algorithm, for various values of the gauge coupling and the magnetic fiux. (To dis- 
cretize the external magnetic field, the smeared links are multiplied by the U(l) links of equation (2.3).) 
For our zero temperature measurements we simulate 24^ x 32, 32^ x 48 and 40^ x 48 lattices, while for 
the finite temperature runs we have lattice configurations with Nt = 6,8 and 10. Finite volume effects 
are studied on the Nt = 6 ensemble using sets of Ns = 16, 24 and 32 lattices. The masses of the up, 
down and strange quarks are set to their physical values along the line of constant physics (LCP) by 
fixing the ratios fK/n^-K and fK/'nn-K to their experimental values. The lattice spacing is determined 
by fx- Details of the determination of the LCP and the lattice scale can be found in, e.g. [54]. Based 
on the reasoning presented in section 4, we use the lattice spacing measurements at T = and B = Q 
to set the scale also at T 7^ and B ^ {). The nonzero value of the magnetic field may of course modify 
e.g. the pion decay constant, just as the temperature can, but this is not important from the aspect of 
matching the lattice quantities at T = 0, -B = to their experimental values (which are also measured 
at T = and B = 0). 

Altogether we generated several hundred to few thousand thermalized trajectories for each (3 
and Ni, (see list of simulation parameters in appendix C), and performed the measurements on every 
fifth one to decrease autocorrelations. The observables presented in section 3 were measured using 
the random estimator method, with 40 random vectors. The production of configurations and the 
measurements were performed on CUDA-capable GPU clusters at the Eotvos University in Budapest 
and on the Bluegene/P at FZ Jiilich. 

We mention here that the staggered formulation of fermions introduces lattice artefacts due to the 
splitting of hadron states into multiplets with different masses [55]. We keep the lowest lying pion state 
at the physical pion mass, while the other members of the pion multiplet are heavier. In the continuum 
limit this mass splitting between the tastes vanishes. However, at finite lattice spacing it can distort 
thermodynamic quantities. To reduce this splitting we apply stout smearing in the fermionic action, 
which is known to significantly reduce taste symmetry violation [54]. 

6 Analysis details 

To study the i?-dependence of the observables of section 3 we scan a wide interval in both the temper- 
ature T and the flux quantum A'^^. The latter is proportional to eB jT'^ , see equation (2.4), so, since the 
transition spreads over a wide temperature region, the physical magnetic field also changes by up to 
a factor of two along an = const, line between T = 120 MeV and T = 180 MeV. To correct for this 
change one can simulate at parameters T, A^f, tuned such that the physical magnetic field remains the 
same. However, since N^, cannot be varied continuously, here we follow a slightly different approach. 
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We measure our observables along a grid of points 
in the T — Nf, plane, as depicted in figure 2. The simu- 
lation points are denoted by the blue crosses, while the 
eB = const, curves are shown by the red dashed lines. 
To perform the interpolation of the measurements along 
these lines in a systematic and effective way, we fit a two- 
dimensional spline function to the data points. A similar 
approach is described in [56] for the fitting of the gradient 
of a two-dimensional function. In figure 3 we show the ob- 
servables as functions of T and for our Nf = 6 lattices. 
We obtain reliable results with good fit qualities; x^/dof. 
being in the range 1.2 — 1.8. 

We perform simulations over the same physical temperature and magnetic field range for two 
smaller lattice spacings at A'^t = 8 and Nt = 10, with very similar x^/dof. values for the spline fits as 
above. We use these three lattice spacings (around Tc(0) they correspond roughly to a = 0.2, 0.15 and 
0.12 fm) to extrapolate our results to the continuum limit. 



1 40 1 60 
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Figure 2: Our simulation points on 24'^ x 6 
lattices (blue crosses) and the lines of constant 
magnetic field (red dashed lines). 




Figure 3: The renormalized up quark condensate (upper left panel), its susceptibility (upper right panel), and 
the strange susceptibility (lower panel) as functions of T and Nf, on our Nt = 6 lattices (note that viewpoints 
are different in order to better show the interesting structures in the particular observables). Measurements are 
denoted by the blue points, while the red surface is the spline fit to the data. The corresponding fit qualities 
are x^/dof. « 1.8, 1.5 and 1.2, respectively. 
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7 Behavior of the condensate 



We remark already at this point that the pseudocritical temperature - as probably best visible in 
the upper right panel of figure 3 for the chiral susceptibility - apparently decreases with increasing 
Nf, ^ B, thereby contradicting a vast number of model calculations, see the summary given in the 
introduction. Furthermore this observation also disagrees with the lattice result of [37]. First of all, 
to check our simulation code we reproduced the results of [37] at a couple of points, see appendix D. 
Since we find a perfect agreement, we conclude that we are left with three possible reasons for the 
discrepancy. First, the lattice spacing of [37] is larger, a ~ 0.3 fm, and also an unimproved action is 
used, so lattice discretization errors may be significant. Second, the present study uses Nf = 2 + 1 
flavors as opposed to the Nf = 2 of [37], and the pseudocritical temperature is known to depend on 
the number of flavors [57], which may also introduce systematic differences in the dependence on the 
external field. Third, the quark masses of [37] are larger than in the present study, which can also 
cause drastic changes in thermodynamics - for example the nature of the transition at B = depends 
very strongly (and non-monotonically) on the quark masses. 
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Figure 4: The unrenormalized chiral condensate as a function of the flux quantum for various temperatures 
around the transition for A'^ = 6. A complex dependence uu{T, TVf,) is observed, since in the deconfined phase 
in some regions the condensate decreases with growing A'^f, (left panel). The magnetic field B'^^^ where the 
renormalized condensate uu^{T, B) is maximal, as a function of the temperature (right panel), as measured on 
Nt = 6 lattices. 



On closer inspection, the differences between our results and those of [37] can actually be traced 
back to the behavior of the chiral condensate as a function of B for a given temperature. While the 
authors of [37] observed that at any temperature the condensate increases with B, we find that this 
dependence is more complex, see the left panel of figure 4 for our Nf = 6 results. At T = 155 MeV, 
which is just above the zero-field pseudocritical temperature, the bare condensate decreases by a factor 
of 2 between A*";, = and Nf^ = 70. As the temperature is reduced the uu{Niy) function starts to develop 
a maximum, clearly visible for T = 142 MeV and T = 136 MeV. This non-monotonic behavior is not 
due to the saturation effects caused by the periodic implementation of the magnetic field on the lattice, 
since this maximum is located at very different values of Nb for temperatures differing only by a few 
percent. Furthermore, for high temperatures the decrease is already visible at < 10 which is in the 
first 5 percent of the period, even for the up quark. To better illustrate this effect and to show that 
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renormalization and conversion from N{y to B does not change the picture quahtatively, in the right 
panel of figure 4 we plot the value of the external field B™^^ where the renormalized chiral condensate 
takes its maximum, as a function of the temperature. At high temperatures this maximum is located 
at B"^^^ = 0, while below T = 155 MeV it shifts to a nonzero magnetic field, in accordance with the 
left panel of figure 4. As already mentioned in the introduction, the possibility of such a decrease in 
the condensate with B was also raised in low energy model calculations [36]. 
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Figure 5: The bare chiral susceptibility for large and vanishing magnetic fields (upper panels), and the dif- 
ference between the up condensates at V eB = 0.9 GeV and at i? = (lower panels) for the iVt = 6 lattices. 
Results are shown for the Nf = 2 + 1 theory (left panels) and for the Nf = 3 theory, where each quark has the 
physical strange quark mass (right panels). 

We summarize our findings as a) the dependence of the condensate on the external field is non- 
monotonic and varies strongly with temperature, and b) as a result the pseudocritical temperature 
shifts to lower values at large B as compared to the B = case. The latter observation is supported 
by a similar Tc{B) dependence deduced from the chiral susceptibility or the strange quark number 
susceptibility, see section 9. We investigate the reason for this behavior further by increasing our 
light quark masses up to the physical strange quark mass, studying the Nf = 3 theory. As a first 
approximation we apply the same lattice scale and line of constant physics as was used for the Nf = 
2 + 1 flavor analysis. This clearly introduces a systematic error, but most probably does not affect 
the qualitative behavior. In figure 5 we show the Nf = 2+1 results for the chiral condensate and 
susceptibility (left panels), compared to the Nf = 3 data (right panels). In the upper panels we plot 
the unrenormalized chiral susceptibility for the case of a vanishing external field (red bands) and a 
large field of ^/ eB = 0.9 GeV (blue bands). Furthermore, in order to see the change in the condensate 
due to the presence of the external field, we plot the difference between the condensate at \/ eB = 0.9 
GeV and at S = in the lower panels. For the case of Nf = 2 + 1 we plot the B = const, slice of the 
2-dimensional surfaces we obtained as described in section 6, while for Nf = 3 we fit the data to a 
simple spline function (in the latter analysis we keep the physical value of the magnetic field constant 
by tuning A^;, ~ B /T"^ as a function of T to keep B fixed). 

As is clearly visible in the lower left panel of figure 5, the magnetic field reduces the chiral con- 
densate for temperatures T > 140 MeV, thus pushing the infiection point of the condensate towards 
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the left and causing a decrease in Tc{B). This decrease is also visible from the behavior of the cor- 
responding susceptibility, shown in the upper left panel of the figure. On the other hand, for larger 
quark masses the situation drastically changes: the condensate increases with the magnetic field for all 
temperatures (see the lower right panel of figure 5), similarly as was observed in [37]. Moreover, there 
is no clear change in T^. the chiral susceptibility (upper right panel) gives consistent pseudocritical 
temperatures for both i? = and ^J eB = 0.9 GeV. This observation supports our explanation number 
three, namely that the difference regarding the change in Tc{B) between the present work and the 
study of [37] stems at least partially from the larger-than-physical quark masses of the latter.^ 

To illustrate the behavior of the condensate from yet 
another aspect, we show in figure 6 the contour plot of 
the renormalized chiral condensate as a function of T and 
B. The color of the curves encodes the value of the con- 
densate, ranging from —0.18 to 0.18 (blue towards red) in 
steps of 0.03. A similar plot of the results of [37] would 
consist of curves having positive slopes, indicating that 
each point of uu moves towards the right as a result of a 
finite B. Here we find that for example the uu^' = —0.12 
curve clearly has segments with negative slope, which once 
again refiects the complex behavior of the condensate as 
a function of B and T. 




130 140 150 
T (MeV) 

Figure 6: Contour plot of uu^ . 



8 Nature of the transition - finite size effects 

Here we address the question of how the strength of the transition changes as the external field is 
switched on. At B = the transition is known to be a broad crossover [58], where the approximate order 
parameters like the chiral condensate change smoothly with the temperature, and no finite volume 
scaling is visible in the observables. Furthermore, the crossover nature of the transition implies that - as 
we will also observe, see figure 9 - different observables give different pseudocritical temperatures [59]. 

As can be seen from figure 3, not just the transition temperature changes with B, but also the 
shapes of our observables as functions of T are altered by a finite magnetic field. Specifically, we 
find that the maximum value of the chiral susceptibility Xu increases with B, which may suggest 
the transition to become stronger for large magnetic fields, as was also reported in [37]. To properly 
determine the nature of the transition we search for finite volume scaling in our observables. 

To this end we perform simulations at our largest magnetic field on the Nt = 6 lattices with 
Ng = 16,24 and 32. The largest lattice in the transition region corresponds to a box of linear size 
~ 7 fm. Here we keep eB/T"^ fixed (and not B itself) as we are only interested in differences between 
the various volumes. In figure 7 the results for the chiral susceptibility (left panel) and for the chiral 
condensate (right panel) are shown as functions of the temperature for eB /T"^ ~ 82. The figure shows 
that our Ns = 16 results agree within statistical errors with the Ng = 24 and Ng = 32 data, indicating 

^Note that while the mass of the Goldstone pion of [37] is below 200 MeV, due to the larger taste splitting the higher 
lying pion tastes may have a larger impact on the response to the magnetic field. 
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Figure 7: The unrenormalized chiral susceptibility (left panel) and chiral condensate (right panel) as functions 
of T measured on our Nt = 6 lattices for different spatial volumes. No finite size effects are visible within 
statistical errors. The crossover nature of the transition persists up to this large external field. 
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that finite size errors are small, compared to statistical errors. This observation also implies that the 
transition at this high magnetic field is still an analytic crossover. 

To further study how the strength of the tran- 
sition changes we investigate the width of the renor- 
malized chiral susceptibility. In figure 8 we plot the 
susceptibility divided by its maximum value as a 
function of T — Tc{B) for three different values of 
the magnetic field for the Nt = 6 lattices. We find 
that, although the height of the peak in Xu gi'ows 
significantly (by almost a factor of 2 between B = 
and the largest B, see also upper left panel of fig- 
ure 5), the width of the peak is only mildly affected 
by the magnetic field. In particular, the width of the pjg^^^ g. p^^j^^j^g ^^^^^^^ T-dependence of 

peak at half maximum decreases from ~ 30(3) MeV Xu as measured on A^t = 6 lattices. The width of the 
to ~ 25(3) MeV as the external field is increased P o Y Y 

from zero to eB = 1.05 GeV^ We find a very similar behavior on the iVj = 8 and 10 lattices. From 
this analysis our final conclusions are that the width of the transition decreases only mildly with 
increasing magnetic field, and as the finite size scaling analysis has shown, the transition remains an 
analytic crossover at least up to V eB ~ 1 GeV. 



- 



1 1 1 1 1 1 

_ M eB=1 .05 


1 1 1 1 1 1 1 


1 1 1 1 


. \K eB = 0.53 


GeV2 #^ 




. m eB = 






~, , 1 , , , 


1 . , , 1 , . 


, 1 , ,~ 



-40 



-20 
T-T^(B) (MeV) 



20 



9 The phase diagram 

Finally, using the fitted two-dimensional surfaces of section 6, we study the observables as functions of 
the temperature, along the lines of constant magnetic field. In particular we analyze the renormalized 
chiral susceptibility X« + X^i the renormalized chiral condensate uu^' + dcT and the strange quark 
number susceptibility c|. For the latter two observables we determine the pseudocritical temperature 
Tc{B) as the inflection points of the curves, while for the former we calculate the position of the 
maximum value of the observable. The results are shown in figure 9. 
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Figure 9: The phase diagram of QCD m the B — T plane, determined from the renormahzed chiral condensate 
uu^ + d(F (upper left panel) , the renormahzed chiral susceptibility Xu + Xd (upper right) and the strange quark 
number susceptibility (lower panel). 

To carry out the continuum extrapolation, we fit the results for Tc{B) for all three lattice spacings 
{Nt =6,8 and 10) together with an A^^-dependent polynomial function of order four of the form 
Tc{B,Nt) = Yli=oi(^i + biNf'^)B\ This ensures the scaling of the final results with iV^"^ ~ a^. We 
obtain x^/dof. ~ 0.5 . . . 1.2 indicating good fit qualities. In order not to make the plots overcrowded, 
we only show error bars for the continuum curves. The error coming from the continuum extrapolation 
is estimated to be 2 MeV and is added to the statistical error in quadrature. The error in the lattice 
scale determination [54] propagates in the Tc{B) function and amounts to an additional 2 — 3 MeV 
systematic error, which is not added to the errors for figure 9 since we find that it does not infiuence 
the shape of the curves. 

As is clearly visible in figure 9, all three observables show that the pseudocritical temperature 
decreases with growing external field B. Preliminary results for the Polyakov loop at one lattice spacing 
show a very similar decrease in Tc{B), see appendix E. We observe that the strange susceptibility (which 
can be viewed as a quantity signaling the deconfinement transition) is the observable most sensitive 
to the external field. Tc{B) changes most drastically in this case, by almost 35 MeV between B = 
and eB ^1 GeV^. We note that our results at i? = are all consistent with earlier determinations of 
the pseudocritical temperature where the stout smeared staggered lattice action was used [51, 59, 60]. 
We mention that one would expect O^a^) effects to become more pronounced as the magnetic field 
grows. However, the numerical data of the pseudocritical temperatures seem to scale well, even up to 
our maximum value of eB ~ 1 GeV^. 
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10 Summary 



In this paper we studied the finite temperature transition of QCD in the presence of external (elec- 
tro)magnetic fields via lattice simulations at physical quark masses. The extrapolation to the con- 
tinuum limit is carried out, and finite size effects are under control. The results are relevant for the 
description of both the evolution of the early universe and of noncentral heavy ion collisions. 

We obtained the phase diagram of QCD in the B — T plane using three observables in the 
phenomenologically interesting region of < eB < 1 GeV^. Performing a finite volume scaling study 
we found that the transition remains an analytic crossover up to our largest magnetic fields, with the 
transition width decreasing only mildly. This rules out the existence of a critical endpoint in the B — T 
phase diagram below eB = 1 GeV^. Moreover, our results indicate that the transition temperature 
significantly decreases with increasing B. This result contradicts several model calculations present in 
the literature which predict an increase in Tc as B grows (see the summary in section 1). We presented 
indications that the response of Tc to the external field can be traced back to the behavior of the chiral 
condensate as a function of T and B. We showed that this behavior is more complex than is predicted 
by most model calculations (where the condensate increases with B for any temperature), and that it 
depends very strongly on the quark masses. 




eqrjy _ 
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Figure 10: Our final result: the QCD phase diagram in the magnetic field - temperature plane. The colored 
bands represent the pseudocritical temperature as defined from inflection points of the renormalized chiral 
condensate uu^ + dcT (red) and the strange quark number susceptibility c| (blue) in the continuum limit. Also 
indicated by the dashed vertical lines are the maximal magnetic fields produced at RHIC and at the LHC. The 
large B region of the phase diagram is relevant for the evolution of the early universe. 

We summarize our results in figure 10, which shows the QCD phase diagram in the B — T 
plane as defined using the renormalized chiral condensate uu^ + dd^ and the strange quark number 
susceptibility c| in the continuum limit. By comparing our magnetic fields to the maximal fields that 
may be produced in noncentral heavy ion collisions we conclude that the decrease in Tc is negligible 
for RHIC and may be up to 5 — 10 MeV for the LHC. Moreover, the effect grows with the magnetic 
field, exceeding 20% for c| aX eB = I GeV^ This may have a significant impact on the description of 
the QCD transition during the evolution of the early universe. 
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A Lattice vector potential and periodic boundary conditions 

In this appendix we show that the lattice prescription for the U(l) links, as in equation (2.3), is indeed 
equivalent to the continuum vector potential up to a local U(l) gauge transformation. The direct 
lattice discretized version of the continuum vector potential (2.1) can be written as 

Uuin) = 1, (z^ / y), 

(A.l 

The periodic boundary conditions are now only satisfied up to a local U(l) gauge transformation 
(transition function), 

Uy{N^,ny,nz,nt) = Uy{0,ny,nz,nt) -V, 

where u is the abelian gauge field and V the gauge transformation that acts in U(l) space. However, 
on the lattice it is more convenient to have exactly periodic boundary conditions. Hence we perform 
the inverse U(l) gauge transformation on the last x-slice of the lattice, which changes fermions ip as 

il){Nx,ny,n:„nt) '4!{N^,ny,nz,nt) ■ , (A.3) 

and the links in both the x and y directions as 

Uy{Nj;,ny,n;,,nt) ^ Uy{N^,ny,nz,nt) -V'^, 

(A.4) 

Ux{Nx - l,ny,nz,nt) u^iN^ - l,ny,n;,,nt) ■ V "'^ 
resulting in periodic boundary conditions and the "twisted" links that we presented in equation (2.3). 



B Renormalization properties of eB from U(l) gauge invariance 

QCD with an external magnetic field has a local U(l) gauge invariance. Let V'o be the bare quark field 
and ^0 a-iid eo the bare external electromagnetic field and electromagnetic coupling. The renormaliza- 
tion of these are given as. 



(B.l) 



If both the regularization and the renormahzation prescriptions are gauge invariant, then so are the 
renormahzation constants Z^^Zj, and Zf.. The gauge transformation for the bare and renormahzed 
quark fields is of the form 

iP' = ^° exp (ia) , ip^' = i)^ exp {ia^) , (B.2) 

which shows that a = , due to the gauge invariance of Z2 . The same transformations for the external 
electromagnetic field are then 

Al' = Al + ^d^a, = Aj^ + ^d^a. (B.3) 

Dividing the second equation by y/Z^ and equating it to the first we obtain (inserting = Z^ ■ e^), 

which implies 

Z,y%=l. (B.5) 

This is a well known result in QED which therefore also applies for the case of QCD with an external 
magnetic field. Since Z(.-\fZz is the particular combination which renormalizes the product eS, our 
conclusion is that eB does not need renormahzation. 



C Simulation parameters 

In this appendix we tabulate the simulation parameters for the T = and the T > runs. The number 
of thermalized trajectories generated for each set of parameters (13, N^) ranges from several hundred 
to a few thousand. 



16^ X 6 


13 


3.45, 3.465, 3.48, 3.488, 3.492, 3.495, 3.497, 3.5, 
3.505, 3.507, 3.51, 3.514, 3.518, 3.525, 3.54 


Nb 


31 


X 6 


f3 


3.45, 3.465, 3.48, 3.495, 3.51, 3.525, 
3.54, 3.555, 3.57, 3.585, 3.6, 3.625 


Nb 


0, 1, 2, 3, 4, 8, 16, 24, 32, 50, 70 


323 y, g 


(3 


3.465, 3.48, 3.488, 3.492, 3.495, 3.497, 3.5, 3.503, 
3.505, 3.507, 3.51, 3.514, 3.518, 3.525, 3.54 


Nb 


124 


243 X 8 


13 


3.525, 3.55, 3.575, 3.6, 3.625, 3.64, 
3.65, 3.675, 3.7, 3.725, 3.75, 3.775 


Nb 


0, 1, 2, 3, 4, 5, 7, 9, 13, 18, 29, 40 


28^ X 10 


13 


3.6, 3.625, 3.65, 3.675, 3.687, 3.7, 3.712, 3.725, 
3.738, 3.75, 3.762, 3.775, 3.8, 3.825, 3.85, 3.875 


Nb 


0, 1, 2, 4, 8, 12, 16, 25, 34, 40 



Table 1: Simulation points for the T > runs. 
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24^ X 32 


P 


3.45, 3.55 




1, 2, 3, 4, 8, 12, 16, 20, 24, 32, 70 


32^ X 48 


P 


3.67 


Nb 


1, 2, 3, 4, 5, 6, 12, 24 


40^ X 48 


/3 


3.75 


Nb 


1, 2, 3, 4, 6, 12 



Table 2: Simulation points for the T = runs. 



D Code check 

To check our code and simulation algorithm we 
reproduced the results of [37] at one tempera- 
ture. We employ exactly the same simulation 
setup, i.e. we use the Wilson gauge action and 
Nf = 2 flavors of unsmeared naive staggered 
quarks on a 16'^ x 4 lattice. We measure the up 
and down quark condensates and the Polyakov 
loop (see definition in appendix E) at gauge 
coupling /3 = 5.35 and mass am = 0.075. We 
plot the average condensate and the Polyakov 
loop in figure 11 (to conform to the notation 

of [37] we divide our Polyakov loops by 3). Figure 11: The average condensate and the Polyakov 

TIT 1 ii i li r 1 ii - , T r 1 loop as functions of the maffnetic flux in the Nf = 2 the- 

We observe that results for both uu + ad and , , r 

ory, compared to results or [37]. 

P agree within statistical errors, with the ex- 
ception of one point for P where the values differ by 2a, as expected for 8 points on statistical 
grounds. Therefore we confirm that there is no discrepancy between results from the two algo- 
rithms/implementations. 

E Polyakov loop 

We carry out the same analysis as presented in section 6 for the Polyakov loop, 

^ n ^4(n). (E.l) 

n^,ny,n, nt=0 

We note that while the quark condensates and susceptibilities and the quark number susceptibility 
depend explicitly on the magnetic field, the Polyakov loop, as a purely gluonic operator, is only affected 
by the modified spatial links indirectly; its expectation value at S > is infiuenced by the magnetic 
factors of equation (2.3) appearing in the fermion determinant. To cancel the multiplicative divergences 
of P, we define the renormalized Polyakov loop [59] using the static quark-antiquark potential V{r) 



Q- 0.2 h 

^ 0.15 
l-o 

-I- 

p 0.1 



-| 1 1 1 1 1 1 1 1 1 1 1 1 r 

• (uu+dd)/2 N,=2 ^ 
(uu+dd)/2 [1005.5365] 



■ P/3 N, = 2 

• P/3 [1005.5365] ! 
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as 



P'^iB, T) = P{B, T)e^(^o)/2T_ 2) 

We measure this observable and perform the 
sphne fitting, see figure 12 for the Nt = G 
results. The inflection point moves to smaller 
temperatures as we increase the magnetic field. 
This behaviour is similar to the decreasing 
transition temperature observed for the con- 
densate or susceptibility (c.f. figure 3). 
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